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The goal of high-energy nucleus-nucleus collisions is to produce and characterize a 
state of nuclear (QCD) matter at (energy) densities well above the nuclear ground state 
(£o ~0.15 GeV/fm 3 ). At high densities and/or at high temperatures one expects [1, 2] 
that quarks are no longer confined in protons and neutrons but move freely over distances 
larger than the size of the nucleon (~1 fm=10~ 15 m). Such a deconfined state of matter, 
earlier named the Quark-Gluon Plasma (QGP) [3], was the state of the Universe within 
the first (~10) microseconds of its creation in the Big Bang [4] and may exist as well 
in the core of neutron stars. The characterization of quark-gluon matter in terms of its 
equation of state (EoS, relating pressure to energy) and of its transport properties (like 
viscosity) and delineating its phase diagram [5] is a major ongoing research effort. See 
[6] for a recent extended treatment of the topic. 

At low energies (beam energies per nucleon of up to 10 GeV/A on fixed target, 
corresponding to center of mass energies per nucleon pair, ^/^nn ^5 GeV) it is expected 
that compressed nucleonic matter is produced. The EoS of nuclear matter [7] at densities 
a few times normal nuclear density (po = 0.17 fm 3 = 2.7 • 10 14 g/cm 3 ), expressed as 
the nuclear compressibility, has relevance for the maximum mass of neutron stars (see 
ref. [8] for a recent overview). 

Employing Quantum Chromo-Dynamics calculations on lattice, a deconfmement 
phase transition for an energy density of about 1 GeV/fm 3 was predicted (see [9] for 
an early review). It was shown [10] that the phase transition at zero baryochemical po- 
tential, is of crossover type, namely with a continuous, smooth, increase of ther- 
modynamic quantities. The value of the (pseudo-)critical temperature, T c , at vanishing 
baryochemical potential (/^) is currently estimated to be in the range 155-160 MeV 
[11, 12]. The existence of a critical point (denoting the end of the first order phase tran- 
sition line, a point where the phase transition is of a second order) is a fundamental 
question, addressed both experimentally [13] and theoretically [14]. 

A nucleus-nucleus collision is a highly dynamical event. One can identify, schemati- 
cally, the following stages of the system ("fireball"): i) initial collisions, occuring during 
the passing time of the nuclei (t pass = 2R/y cm c); ii) thermalization: equilibrium is estab- 



lished; iii) expansion and cooling (in a deconfined state); iv) chemical freeze-out (pos- 
sibly at hadronization): inelastic collisions cease, hadron yields (and distribution over 
species) are frozen; v) kinetic freeze-out: elastic collisions cease, spectra and correla- 
tions are frozen. 

The challenge is to characterize the hot (deconfined) stage iii), while most of the 
measurements are performed via hadrons (or their decay products) carrying information 
from the system at stages iv) and v). Even though the early stage of hot deconfined mat- 
ter remains inaccessible in a direct way because of quark confinement, there are experi- 
mental observables which carry information from this stage. Extracting the properties of 
the deconfined stage is possible only via models. At lower energies hadronic transport 
models [7] are employed, while at higher energies hydrodynamics [15] is widely used. 
Based on model comparison to data, one can extract the following ranges of the fireball 
characteristics: Temperature: T = 100 — 1000 MeV, or up to a million times the temper- 
ature at the center of the Sun (1 MeV~10 10 K); Pressure: P = 100 - 300 MeV/fm 3 (1 
MeV/fm 3 ~ 10 33 Pa); Density: p = 1 — 10 • po; Volume: several thousands fm 3 ; Dura- 
tion: 10-20 fm/c (or about 3 - 6 • 10~ 23 s). 

The experimental "control parameters" are: a) the collision energy (per nucleon pair, 
v /«nn); b) the centrality (impact parameter, b) of the collision (or, alternatively, the size 
of the colliding nuclei), which is deduced from specific measurements and involves 
calculations within the Glauber model [16]. A usual way of expressing centrality is 
via the number of participating nucleons, N part , namely the nucleons involved in the 
creation of the fireball in the overlap region of the two colliding nuclei. 

After the initial measurements at the Bevalac (Berkeley) in the '80s, the program 
of heavy-ion collisions continued at higher energies at Brookhaven at the Alternating 
Gradient Synchrotron (AGS) and at CERN at the Super-Proton Synchrotron (SPS), 
while in the low energy range measurements were performed at GSI Darmstadt at 
the Schwerionensynchrotron (SIS). Started in year 2000, the experimental program at 
the Relativistic Heavy Ion Collider (RHIC) at Brookhaven spans ^/snn from ~8 to 
200 GeV (see, for earlier results, experimental summaries in [17, 18, 19, 20] and an 
overview in [21]). The study of QCD matter has entered a new era in year 2010 with the 
advent of Pb-Pb collisions at the Large Hadron Collider (LHC), delivering the largest 
ever collision energy, v /snn=2.76 TeV, more than a factor of 10 larger than previously 
available. A recent overview of the first LHC data is available in ref. [22]. 

The measurement of the charged hadrons pseudo-rapidity (r] = — ln[tan(0/2)], with 
the polar emmission angle) density, dA^h/dr/, at the LHC [24] was eagerly awaited 
and showed that the increase compared to the measurement at RHIC is by a factor of 
about 2.2 for central collisions. Interpreted as the outcome of an increase of the initial 
entropy density, this increase can be translated into a factor of 1.3 increase of the initial 
temperature [25]. The ALICE measurement confirmed the phenomenological (^/snn) 0,3 
behavior seen at lower energies [26], see Fig. 1 (left panel). Shown are rapidity densities 
of charged particles dN c h/dy at mid-rapidity, y=0 (where particles are emitted in the 
transverse direction); y = \ In jtzj^ = tanh^pY), with (pY) the longitudinal (beam 

direction) momentum (velocity), E = ^/m 2 + p 2 the total energy. The measurement at 
the LHC clearly demonstrated that the increase of dA^h/dr] with energy is steeper in 
nucleus-nucleus (AA) collisions compared to pp collisions, where the functional form 




FIGURE 1. Left panel: collision energy dependence of charged particle density cW c h / dy at midrapidity, 
measured by various experiments in central collisions corresponding to Af parf =350. Right panel: energy 
dependence of the transverse energy (plot from [23]). 



is (^/snn) [24]. The centrality dependence of dA^h/dr/ is at the LHC [27] identical 
to that measured at RHIC, pointing out to a similar mechanism of particle production 
at the two energies. A model of the parton structure of matter at low parton fractional 
momentum x, the Color Glass Condensate [28], is in a good agreement with the data 
(see an extended comparison to theoretical models in [27]). The data points shown in 
Fig. 1 (left panel) are obtained by summing the measured dN/dy yields for pions, kaons 
and protons and their antiparticles, see below. 

Utilizing, in addition to particle counting, the momentum measurement (or, alterna- 
tively measuring the total hadron energy in calorimeters), see Fig. 1 (right), one can ex- 
tract from the data the energy density at the thermalization time. This involves a space- 
time model of the collision, which was put forward by Bjorken [29]. In this model, 
the energy density is e = = ^^H"' wnere is the transverse energy and 

At = TtR 2 is the geometric transverse area of the fireball (for central Pb-Pb collisions, 
Aj ~150 fm 2 ). Assuming a conservative value for the equilibration time, Tq=\ fm/c, one 
calculates for the LHC energy a matter energy density of £l#c=15 GeV/fm 3 [23], well 
above the "threshold" value for quark-gluon matter, calculated in lattice QCD [9], of 
about 1 GeV/fm 3 . 

In Fig. 2 (left panel) the collision energy dependence of identified hadron yields at 
mid-rapidity is shown. This comprises measurements by experiments at the AGS: E895 
[30, 31, 32], E866/E917 [33, 34], E891 [35]; the SPS: NA49 [36, 37, 38, 39], NA44 [40], 
NA57 [41]; RHIC: STAR [42, 43, 44, 45], BRAHMS [46], PHENIX [47]; LHC: ALICE 
[48]. The monotonic decrease of the proton yield as a function of energy indicates that 
fewer and fewer of the nucleons (or their valence u, d quarks) in the colliding nuclei are 
"stopped" in the fireball. An onset of meson production is seen, with the kaons (heav- 
ier and containing a strange quark) produced less abundantly than pions. The asymmetry 
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FIGURE 2. Left panel: hadron multiplicities at mid-rapidity in central collisions. Right panel: the phase 
diagram of strongly interacting matter with the points representing the thermal fits of hadron yields. The 
down-pointing triangle indicates ground state nuclear matter (atomic nuclei). 



between the n + and n~ production yields reflects the isospin composition of the fireball. 
The asymmetry between K + and K meson and A and A hyperon production is deter- 
mined by the quark content of the hadrons: K + (us), K~{us) A(uds), Muds). At lower 
energies, the availability in the fireball of valence u, d quarks from colliding nucleons 
leads to a preferential production of hadrons carrying those quarks. These asymmetries 
vanish gradually for higher energies, where the hadrons are mostly newly created (re- 
flecting Einstein's famous equation m = E/c 2 ) and the production yields exhibit a clear 
mass ordering. Good fits of the measurements are achieved with the thermal model with 
3 parameters: temperature T, baryochemical potential \x b , and volume V. The thermal 
model describes a snapshot of the collision, namely the chemical freeze-out, which is 
assumed to be quasi-instantaneous. This provides a phenomenological link of data to the 
QCD phase diagram, shown in Fig. 2. Each point corresponds to a fit of hadron yields in 
central collisions of Au-Au or Pb-Pb nuclei at a given collision energy. A remarkable 
outcome of these fits is that T increases with increasing energy (decreasing /if,) from 
about 50 MeV to about 160 MeV, where it exhibits a saturation for fi b <300 MeV. This 
saturation of T led to the connection to the QCD phase boundary, via the conjecture that 
the chemical freeze-out temperature can be the hadronization temperature [49, 50]. This 
picture is currently tested at the LHC with first identified hadron yields data [48]. 

Collective flow is a distinct feature of nucleus-nucleus collisions. In central collisions 
one investigates the so-called radial flow, which is quantified fitting transverse momen- 
tum (pt) spectra with the so-called "blast wave" model [51], accessing in a convenient 
(albeit maybe over-simplified) way bulk properties of the fireball at kinetic freeze-out. 
The extracted fit parameters, the temperature T and average transverse velocity (j8r)» 
are shown in Fig. 3 (left) as a function of the collision energy. The measurements are 
by experiments FOPI [52], EOS [53], NA49 [54], STAR [42], ALICE [48]. A strong 
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FIGURE 3. Collision energy dependence of collective flow: Left panel: radial flow in central collisions, 
quantified by the temperature T and average velocity (j5j) at kinetic freeze-out; the temperature at 
chemical freeze-out is represented by the triangles. Right panel: elliptic flow in mid-central collisions. 



increase of both T and (j8t) is seen at low energies (beam energies of up to 1 GeV/A 
on fixed target) with a small further increase of (jSp) and a constant kinetic freeze-out 
T, which is 50-60 MeV lower than the chemical freeze-out T. At lower energies, the 
chemical freeze-out T is smaller than the kinetic T, which is unphysical and awaits a 
resolution. At the LHC, (jSp) ~0.65c [48]. 

In Fig. 3 (right) we show the energy dependence of elliptic flow, measured in mid- 
central collisions (N part — 160, corresponding to an average impact parameter value 
of (b) ~ 7 fm) by experiments FOPI [55], E895 [56] E877 [57], STAR [58], ALICE 
[59]. Elliptic flow arises in non-central collisions of nuclei as a result of the initial 
elliptic trasverse shape of the overlap zone of the two nuclei (participant eccentricity 
£part)- Through the initial gradients of the energy density (or pressure), this leads to 
anisotropic spatial (angular) emission of hadrons. This is quantified by the second order 
(quadrupole) Fourier coefficient v% = (cos(20)), where is the azimuthal angle with 
respect to the reaction plane. The complex evolution of elliptic flow as a function of 
energy seen in Fig. 3 (right) is understood qualitatively rather well. At low energies, in- 
plane (v2 > 0), rotation-like, emission may arise due to low energy density in the overlap 
region and of long reaction times. The fast transition towards preferential emission 
out-of-plane (v2 < 0) is the outcome of more energetic collisions, leading to a larger 
energy density of the fireball. The increase of elliptic flow is a fingerprint of a stronger 
collective expansion, hindered by the passing spectators, which act as a shadow for the 



outgoing nucleons and fragments. The competition between the increasing speed of the 
expansion and of the decreasing passage time t pass of spectators leads to a maximum 
of (absolute value) elliptic flow in the SIS energy range. In this energy domain, the 
transiting spectators, with t pass varying between 40 and 10 fm/c, act as a clock for the 
collective expansion. In this regime, elliptic flow (historically called "squeeze-out") is a 
prominent observable for the extraction of the nuclear EoS [7]. Towards larger energies, 
elliptic flow exhibits another transition, to a preferential in-plane emission [60], a result 
of a unhindered collective expansion of the initially-anisotropic fireball. Elliptic flow is 
built mostly in the earlier stages of the collision, since it is determined by the initial 
pressure gradients, which it alters quickly as it develops. Consequently, at high energies, 
elliptic flow probes (albeit not exclusively) the deconfined state of quarks and gluons. 

The long-awaited elliptic flow measurement at the LHC [59] exhibits a larger mag- 
nitude by ~35% compared to the measurement at ^/snn = 200 GeV. This increase 
is described by hydrodynamics and was anticipated on a purely phenomenological 
log( v /5NN) behavior seen at lower energies [26]. The data show [59] that the increase 
is due exclusively to the increase of the average transverse momentum of the hadrons, 
while the pj dependence of V2 is identical at LHC to that measured at RHIC. 

Another set of experimental observables from the kinetic freeze-out stage are obtained 
from Hanbury Brown and Twiss (HBT) interferometry, which allows to measure to the 
size of (a region of) the fireball and its lifetime [61]. Elliptic flow, radial flow and HBT 
interferometry are the observables used to extract, via comparisons to hydrodynamic 
calculations, the viscosity of deconfined matter (more precisely, shear viscosity divided 
by entropy density, r\/s [62]). The remarkable description of flow and HBT interferom- 
etry in hydrodynamic models observed at RHIC is further confirmed and extended with 
the measurements at the LHC. The quantitative determination of rj/s is dependent on 
the initial energy density distributions, which are calculated according to either Glauber 
(binary collisions scaling) or Color Glass Condensate model. It appears that deconfined 
matter is characterized by low values of r\ fs; for an estimate for nuclear matter in mildly- 
excited nuclei see [63]. Thermalization at a very early stage of the collision, at or below 
1 fm/c, as needed in hydrodynamic calculations, is presently a challenge to theory [64]. 

Can we identify an onset of deconfinement based on (bulk) hadronic observables 
discussed above? The change, in the range ^/snn —5-10 GeV, of fireball properties 
(see Fig. 3) is a possible fingerprint, but further experimental and theoretical support 
is needed to conclude. 

Probing the deconfined matter in a more direct way is done with special observables 
of the early stage. One category comprises thermal photons and low-mass dileptons [65], 
produced in the hot (deconfined) matter over its entire lifetime. Measurements of thermal 
photons at RHIC have shown [66] that the temperature averaged over the lifetime of the 
fireball is larger than the chemical freeze-out T . Another category of QGP probes is the 
so-called hard probes, namely processes characterized by an energy scale (quantified by 

the transverse mass mj = \Jm^-\- p\, where mo is the rest mass of the hadron) above 
several GeV (well above the temperature of the medium). Examples of such observables 
are hadrons at high pj (or jets) and hadrons containing heavy (charm or bottom) quarks. 
They are produced at early times in the collision (t = l/mj, with 1 fm = 1/(0.2 GeV) in 
the system of units where h = c = 1 commonly used in high-energy physics). 



Proposed by Bjorken in 1982 [67], "jet quenching", the extinction of jets (due to 
the energy loss of the parent parton) in QGP was for the first time observed at RHIC 
[17, 18, 19, 20] and is a subject of intense study at the LHC [68, 69, 70, 71]. The 
usual method to quantify jet quenching is via the nuclear modification factor, defined 
as: Raa = (^AA/dy) ■ d/V pp /dy), where dN/dy denotes the yield of a given 
observable measured in nucleus-nucleus (AA) or proton-proton (pp) collisions and N co u 
is the average number of nucleon-nucleon collisions over the given centrality interval of 
AA collisions; N co u is calculated using the Glauber model [16]. 

A change of physics in AA collisions (which in specialized terms is called a "medium 
effect") is seen as a departure of Raa from unity. However, modifications of parton dis- 
tributions in nuclei compared to pp (shadowing or saturation) need to be considered 
carefully, in particular at LHC energies; very recent measurements in p-Pb collisions 
address this issue [72]. We note that the binary collision scaling assumed in the con- 
struction of Raa only applies to hard processes. It is known experimentally that bulk 
particle production (comprising essentially pions, protons and kaons at low-momentum, 
Pt ^3-4 GeV/c) in AA collisions scales (in first order) with N part [27]. 
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FIGURE 4. Left panel: the nuclear modification factor 7?aa as a function of transverse momentum 
for charged hadrons in Pb-Pb collsions at the LHC (Fig. from [70], see refs. therein for the theoretical 
curves). Right panel: centrality dependence of Raa for charged hadrons at high pj in comparison to that to 
charmed mesons (D), charmonium (J/y/) and beauty hadrons (B). The dotted lines denote the contribution 
of production in the corona of colliding nuclei. 

The first measurement of Raa for charged hadrons at the LHC [68] showed that the 
suppression is larger than previously measured at RHIC, reaching a factor of about 
7. The suppression is reduced towards larger pj values, see Fig. 4 (left), but remains 
substantial even at 50-100 GeV/c [69, 70]; it is also seen in reconstructed jets [71]. 

The basic features seen in the data are reproduced by models implementing parton 
energy loss in deconfined matter. Towards extraction of transport coefficients, the de- 
scription of jet quenching in theoretical models remains a challenging task [73]. The 
measurements, in conjunction with theoretical models, clearly demonstrate that partons 
lose energy in the hot and dense quark-gluon matter, which leads to broader jets [74]. 

For the first time measured at the LHC in a direct way, the nuclear suppression of 
charmed [75] and beauty [76] hadrons at high pi is shown as a function of centrality in 



Fig. 4 (right). The theoretical expectation is that heavy quarks (charm and bottom) lose 
less energy (by gluon radiation) compared to lighter (up, down, strange) ones [77]. This 
expectation is exhibited by the data, see Fig. 4 (right), although a definite conclusion 
needs further support from experiment. The energy loss suffered by energetic heavy 
quarks in QGP is indicative of their "strong coupling" with the medium, dominated 
by light quarks and gluons. The measurements at the LHC consolidate and extend the 
observation at RHIC of heavy quark energy loss and flow [78]. 
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FIGURE 5. Left panel: centrality dependence of the nuclear modification factor for inclusive J/y/" 
productions at RHIC and LHC energies. Data are compared to statistical hadronization model calculations. 
Right panel: transverse momentum dependence of the nuclear modification factor of J/y/ measured at the 
LHC by ALICE (preliminary data) in comparison with transport model calculations (plot from [79]). 

Among the various suggested probes of deconfinement, charmonium (cc) states plays 
a distinctive role. J/y/ is the first hadron for which a clear mechanism of suppression 
(melting) in QGP was proposed early on, based on the color analogue of Debye screen- 
ing [80]. In the statistical hadronization model [81], the charm quarks produced in initial 
hard collisions thermalize in QGP and are "distributed" into hadrons at chemical freeze- 
out. All charmonium states are assumed to be not formed at all in the deconfined state 
but are produced, together with all other hadrons, at chemical freeze-out. See [82] for 
recent predictions of this model. Kinetic recombination of charm and anti-charm quarks 
in QGP [83] is an alternative quarkonium production mechanism. In this model (see 
[84, 85] for recent results), a continuous dissociation and regeneration of charmonium 
takes place in the QGP over its entire lifetime. 

The measurement of J/y/ production in Pb-Pb collisions at the LHC was expected 
to provide a definitive answer on the question of (re)generation. The data measured at 
high-px [86] show a pronounced suppression of J/i// in Pb-Pb compared to pp collisions 
and of the same magnitude as that of open-charm hadrons, see Fig. 4 (right). This may 
indicate that the high-px charm quarks that form either D or J/y/ mesons had the same 
dynamics, possibly a thermalization in QGP and a late hadronization. 

The first LHC measurement of the overall (inclusive in pj) production [87], showed 
Raa values significantly larger than at RHIC energies, see Fig. 5 (left). The data are well 
described by both statistical hadronization model [82], shown in Fig. 5 (left), and by 
transport models [84, 85], which also describe the pi dependence of Raa [79], Fig. 5 



(right). The transport model calculations show that, as expected, that regeneration is pre- 
dominatly a low-px phenomenon. Discriminating the two pictures, statistical production 
at chemical freeze-out and regeneration during the QGP lifetime, may help providing an 
answer to fundamental questions related to the fate of hadrons in a hot medium [88]. 
Recent measurements at the LHC of production of bottonium (bb) states [89] adds one 
more flavor to the QGP research field. 
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